Read system data and assign an initial approx. solution, or use flat-start 



Form Nodal Admittance Matrix, Initialize iteration count r=l 



Initialize node count p=l, and maximum real and imaginary corrections DEMX=0.0 and DFMX=0.0 




yes 




Compute Q p (r+!) by relation (6) after adjusting complex 
voltage for specified value by relation (10) 



1 


no 

f 


Set SHp = PSH p +jQ p (r+1) 




SetSHp^PSHp+jQp 1 



SetSHp =PSH P +jQ p m 



Compute V p (rf !) from relation (5), and adjust complex 
voltage for specified value by relation (10) 



Compute V p (r+1) from relation (5) 



8 



10 



Compute changes in the real and imaginary parts of the node-p 
voltage by using relations (7) and (8), and replace current value of 
DVRMX and DVTMX respectively in case any of them is larger 



11 



Calculate accelerated value of V ? ir * !) by relation 
(9) and update voltage for immediate use in the 
next node-voltage calculation by V p r = Vp (rH) 





p=p+t 








« — 






( STOP ) 



Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 
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Fig. la: Prior Art: flow-chart of Gauss-Seidel Loadflow (GSL) Method 
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a Read system data and assign an initial approx. solution, or use slack-start 



Form Nodal Admittance Matrix & Initialize ITRP=ITRQ=r=0 




yes 



> CosO>p=Cos (-48) & Sin<P p =Sin (-48) 



Form, Factorize, and Store constant gain matrices [Y8] & [YV] for SSDL-YY method 



>■ Compute APi r &AQ r , where i=l,..,(m+k) and j=l,..,m 




yes 



Compute [RP] using eqn. (17) for PQ-nodes and using eqns. (21) and (22) for PV-nodes 

g 



Solve eqn. (13) for [A6] and update voltage angles using, [9] = [6] + [A9] 



Set voltage magnitudes of PV-nodes equal to the specified values, 
and Increment ITRMTRP+1 and r=(iTRP+rrRQ)/2 

■ ■■ 



Compute AP| r & AQ r , where i=l ,..,(m+k) and j=l ,..,m 



Calculate line flows 
and, output indications 
of transformer tap 

position and reactive I 

power generation at 
PV-nodes 




c 



STOP 



yes 



j Compute [RQ] using eqn. (18) for only PQ-nodes 



Solve eqn. (15) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV] 
. While solving (15), skip all the rows and columns corresponding to PV-nodes 



In case of violation of Q-Iiirrit of any PV-node generator and tap-position limit of any 
tap-changing transformer, adiust the loadflow solution by restricting to violated limits 



ml 



Increment ITRP=ITRP+1 and r=(ITRP+ITRQ)/2 



Fig.lblPrior Art: Flow-chart of Super Super Decoupled Loadflow (SSDL) Method 
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Fig. 2a: One-lino- diagram of IEEE 14-node network 




Fig. 2d Non-redundant Level-l sub-networks of fig. 2b are regrouped to reduce 
the number of processors required without increasing the number of 
nodes in any regrouped sub-network larger than the original largest 
sub-network of 6-nodes 
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Fig. 2b: Level-l sub-networks around circled nodes for the network 
of fig. 2a 
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C START J 



1 Read system data and assign an initial approx. solution, or use flat-start 



Form Nodal Admittance Matrix, Initialize iteration count r=l 



3 Initialize node count p=l , and maximum real and imaginary corrections DEMX=0.0 and DFMX=0.0 4 




yes 



yes 



Compute Qp (r+!) by relation (6) after adjusting complex 
voltage for specified value by relation (10) 



SetSHp = PSHp+jQ p (rH) 

i. 




SetSHp = PSHp + jQp' 



SetSHp=PSH p +jQp m 



Compute V p (rfl) from eqns. (27), (28), and (29) and adjust 
complex voltage for specified value by relation (10) 



Compute V p (r+,) from eqns. (27), (28), and (29) 



88 



10 



Compute changes in the real and imaginary parts of the node-p 
voltage by using relations (7) and (8), and replace current value of 
DVRMX and DVTMX respectively in case any of them is larger 



99 



11 



Calculate accelerated value of V p (rH) by relation 
(9) and update voltage for immediate use in the 
next node-voltage calculation by V p - V p (rH) 



> r=r+l 



p-p+1 




Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 



14 



Fig. 3a: Invention: flow-chart of Gauss-Seidel-Patel Loadflow (GSPL) Method 
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34 



Read system data and assign an initial approx. solution, or use flat-start, arid read J in (Jf, Je) iteration scheme 



Form Nodal Admittance Matrix, Initialize iteration count r^ITRF^ITRE^O, DEMX=DFMX=0.0 



13 



Initialize node count p=l , and maximum imaginary part of voltage correction DFMX=O.0 



From step-44 on 
the next page 




Compute Qp (r+,) by relation (31) after adjusting complex 
voltage for specified value by relation (10) 





no 

r 


Set SHp-PSH p +jQp (r+,) 




Compute f , <If "from eqns. (44) or (46) and (29) and adjust 
complex voltage for specified value by relation (10) 



Compute f p (rH, from eqns. (44) or (46), and (29) 



36 



3| 



Compute changes in the imaginary part of the 
node-p voltage by using relations (50), and replace 
current value of DFMX in case it is larger 



To Step-34 on 
the next page 



Calculate accelerated value of fp (r+,) by relation 
(48) and update voltage for immediate use in the 
next node-voltage calculation by f p r = f p (r+I) 



p=p+l 




33 



> ITRF=ITRF+i 




Q STOP 



Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 



Fig. 3ai Invention: flow-chart of Gauss-Seidel-Patel Loadflow (GSPL) Method 
[given for (Jf, Je) iteration scheme-3: ITRF] 
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From step-33 on the 
previous page 



3 1\ Initialize node count p=l , and maximum real part voltage correction DEMX=0.0 < 




Compute Qp (pf,) by relation (31) after adjusting complex 
voltage for specified value by relation (10) 



3<\ , 


no 

r- 


Set SHp = PSH p +jQp (r+1) 




Set SHp = PSHp+jQ p m 



Compute c p (rfl) from eqns. (43) or (45) and (29) and adjust 
complex voltage for specified value by relation (10) 



3<? 



SetSHp=PSH p + jQp n 
I 



Compute e p <r *" from eqns. <43) or (45), and (29) 



ho 



H2- 



Compute changes in the real part of the node-p 
voltage by using relations (49), and replace current 
value of DEMX in case it is larger 



To step-23 on the 
previous page 



Calculate accelerated value of e p (rf l} by relation 
(47) and update voltage for immediate use in the 
next node-voltage calculation by e p r = e P <r+I) 



p=p+l 




ITRE=ITRE+1 




TIC 



Q STOP J 



Calculate line flows, and 

output indications of 
transformer tap position 

and reactive power 
generation at PV-nodes 



Fig. 3b: Invention: flow-chart of Gauss-Seidel-Patel Loadflow (GSPL) Method 
[given for (Jf, Je) iteration scheme-3: ITRE] 
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SERVER-PROCESSOR 



( START J 



A-PARALLEL-PROCESSOR 



Read system data, assign an initial approx. 
solution, or use flat-start to initialize vector 
V(I), and initialize iteration count r=l 




Form Nodal Admittance Matrix for a sub- 
network, and store it in the local memory 
of an assigned processor 


► 



S3 



Initialize maximum real and imaginary 
corrections DEMX=0.0 and DFMX=0.0, and 
W(I>= CMPLEX (0.0, 0.0) 




read relevant elements of V(I) from 
commonly shared memory by all processors 



Perform two GSPL iterations for each sub- 
network. Steps-4 to 88 or 99 in Fig. 3a to 
calculate new estimate of voltage for each 
* sub-network nodes 



'ST 



Calculate [W(I>V(I)], and Compute changes in the real 

and imaginary parts of the node-p voltage by using 
relations (7) and (8), and replace current value of DEMX 
and DFMX respectively in case any of them is larger. 



Calculate accelerated value of V V(I) by relation 
(9) and Perform operation V(I)=W(I) 



r=r+l 




Co 



Calculate line flows, and 

output indications of 
transformer tap p<Jsition 

and reactive power 
generation at W-nodes 



^ STOP J 



Fig. 3 c: Invention: flow-chart of Parallel-Gauss-Seidel-Patel Loadflow 
(PGSPL) Method 
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Main/Server 
Processor, 
which is a 
delegating 
and 

coordinating 
processor 



* ► 



Processor- 1 



< — ► 



Processor- 1 
Memory 



< ► 



Processor-2 



« — ► 



Processor-2 
Memory 



« 



Processor-n 



« — ► 



Processor-n 
Memory 



Processor- Main 
Memory 



Common 

memory 

shared by all 

processors 

including 

Main/Server 

processor 



Fig. 4: Invented Parallel computer Architecture/organization 
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10 


Obtain on-line data 




20 


f 




Initial Established Specification of Controlled Parameters (Real and Reactive Power 


at PQ-nodes, Real Power and Voltage Magnitudes at PV-nodes, and Transformer 




turns ratio etc.) 





30 



40 



Load-Flo\*£. 
Calculation 



Estimation 




of 




Power System State 




50 - 





NO 




60 



Evaluation, Good? 
(No over load & No 
under/over Voltage?) 

70 



YES 




Correction in One or 
More Controlled 
Parameters 



Effect Actual Control 
To set finally obtained 
Controlled Parameters 



Fig. Si Load-Flow Calculation in Security Control of Electrical Power System 
(PRIOR ART) 
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60 



Impedance 




Breaker 




Admittance 


Network 
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Status 
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Matrix 
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Assign Transformation 
Ratio 



30 



Feeder VARS 



40 



42 



Feeder MW 



44 



Generator Loads 



Adjust Node 

MW 
Distribution 



50 



Assign Initial Approx. or Guess 
Loadflow solution, and Read data of 
Operating Limits of various Network 
Components, including Generator 
Excitation and Transformer Tap- 
changing Limits 



Machine 
Excitation Control 

Transformer 
Tap Control 




Form & Store 
Constant Gain 

Matrices, 
[Y0] & [YV] 
which is not 
required for 
GSPL& 
DGSPL 
Methods 



Perform Loadflow 

Calculation by 
Defining/Forming 
& Solving any of 
the SSDL, GSPL, 
or DGSPL Models 
of the Power 
Network 



64 



Load Flow Solution including 
Generators Reactive Power and 
Transformers Turns Ratio 
Indications 



Fig. 6: Load-Flow Calculation for Voltage Control in Electrical Power 
System (Prior Art) 
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r Tt b 



T2- 

■9- 



FH3 



Fig* 7: An Exemplary 6-node Power System 

Nodes: i, 2,3 are PQ-nodes 

Nodes: 4 and 5 are PV-nodes 

Nodes: 6 is the slack/swing/reference node 

Transformers Tl and T2 are tap-changing 
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